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ABSTRACT 

We present a new, high resolution HST/GHRS spectrum of quasar Q0302- 
003, and use the He + Lyman-o; absorption, together with a high resolution Keck 
spectrum of the HI Lyman-a forest, to probe the distribution and ionization 
state of foreground gas just below the quasar redshift z « 3.3. Within « 4000 
km/sec of the quasar redshift the spectrum shows a substantial flux (r « 1) 
with ll He + Lyman-a forest" absorption features correlated in redshift with the 
HI Lyman-a forest; the absorption in this region is accounted for entirely by 
the discrete components of the forest, indeed the main "Gunn- Peterson edge" 
can be identified with a particular complex of HI absorbing clouds. We attribute 
the lack of continuous absorption from diffuse gas to the "proximity effect" in 
this region, a large bubble where helium is highly ionized by the quasar, and 
use its size to estimate the background flux at the He + ionization threshold. 
The near-quasar data also lead to constraints on diffuse gas density near the 
quasar, tied to the observed quasar flux, and helium abundance, tied to the 
observed quasar spectrum. Far from the quasar redshift, the spectrum displays 
He + absorption (r > 1.3) even in several redshift intervals with no detectable 
HI absorption, implying a soft ionizing spectrum as well as absorption from 
gas between detected HI clouds. The smoothed spectrum displays residual 
flux everywhere with an average optical depth tqp < 2, which indicates a 
low density of redshift-space-filling gas; using constraints from the HI ionizing 
spectrum we estimate f2 9 < 0.01(/i/0.7)~ L5 , and infer that the helium is 
already mostly doubly ionized by this epoch. Our estimates are consistent with 
ionization models based on observed quasar populations, previous limits from HI 
Gunn-Peterson studies, and simulations of the gas distribution in CDM models 
of galaxy formation. 
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1. Introduction 

The line of sight to the z = 3.285 quasar Q0302-003 (hereafter, Q0302) provides a 
unique probe of the intergalactic medium at high redshift. It is the highest redshift quasar 
yet discovered where the light is unobscured down to below 304 A in the rest frame, 
allowing measurement of He + Lyman-a absorption by foreground gas. Besides confirming 
the abundant primordial helium predicted by the Big Bang model, the helium absorption 
records new information about the ionization history of intergalactic gas, especially useful 
for disentangling the roles of stars and quasars in reionizing the universe. Most significantly, 
its relatively higher optical depth provides a better tool than HI for measuring absorption 
by diffuse gas which fills the space between the galaxies and between the identified HI 
"Lyman-o; forest" clouds, in the intergalactic and protogalactic media. Resolved helium 
absorption provides a direct "nonlinear map" of the gas distribution in the most rarefied 
gas occupying most of the volume of space. 

Ideas about the cosmic gas distribution have sharpened quantitatively in recent 
years, due to hydrodynamic simulations which accurately predict the motion of matter 
in hierarchical models of galaxy formation (Cen et al. 1994, Hernquist et al. 1995, 
Miralda-Escude et al. 1996, Rauch et al. 1997, Croft et al. 1997, Zhang et al. 1997; see 
also Bi and Davidsen 1997). Departing from earlier models based on isolated clouds with 
symmetric geometries such as spheres and slabs, simulations of the conversion of uniform 
gas into condensations reveal a dynamical system with a complex geometry where the 
distinction between diffuse gas and clouds is blurred and is not always reflected in the 
appearance of an absorption spectrum. Simulated spectra reveal that gas in the most 
underdense regions, filling the bulk of the spatial volume, is so highly ionized that it 
produces absorption features with very low HI Lyman-o; optical depth. The most abundant 
ion He + however produces optical depths of the order of unity even in these regions, so its 
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absorption is easily detectable, mapping the distribution of cosmic baryons at the lowest 
densities. 

Absorption by He + is also the most direct probe of the hard ultraviolet cosmic radiation 
field, which can be predicted from semiempirical models based on observed quasar and 
absorber populations (Haardt and Madau 1996). The spectral shape also influences other 
observables such as the ratio of CIV to SilV (Songaila & Cowie 1996, Giroux & Shull 
1997, Savaglio et al 1997), so information from helium absorption allows information about 
relative C and Si abundances to be derived. In situations where the ionizing spectrum is 
known, such as the near proximity of a quasar, He + absorption can be compared to HI 
absorption to extract independent information about the primordial abundance of helium, 
an important test of Big Bang Nucleosynthesis. 

The first detection of cosmic He + absorption was made in Q0302 by Jakobsen et al. 
using the Hubble Space Telescope Faint Object Camera (FOC). They found an absorption 
edge and a large "Gunn- Peterson" continuous optical depth, r > 1.7, attributed to He + 
Lyman-a absorption by diffuse gas. A similar observation has also been made of the 
z = 3.185 quasar PKS 1935-692, with a similar result (a lower limit on the optical depth 
t > 1.5, also at 90% confidence, by Tytler & Jakobsen 1996; note that this result is based 
on new data which modify earlier conclusions on this object by Tytler et al. 1995). These 
data however had important limitations; the low (10 A) resolution of the FOC could not 
resolve features from the known HI clouds or even from the larger gaps between them, and 
the inaccuracy of the calibration (20 A) could not place the location of the edge accurately 
relative to either the quasar or the clouds. 

A significant improvement came from the Hopkins Ultraviolet Telescope (Davidsen et 
al. 1996), which can reach shorter wavelengths and hence lower redshift than HST, and also 
provides better resolution and wavelength calibration than FOC. Davidsen et al. observed 
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the z = 2.72 quasar HS 1700+64 and found an He + edge close enough to the predicted 
redshift to rule out the possibility of foreground HI as an important contaminant. They 
also found that the flux below the edge is not consistent with zero, and measured accurately 
a mean optical depth, r = 1.00 ± 0.07. The decreasing absorption with time reflects the 
increasing ionization of He + at around z = 3 and the conversion of diffuse gas into clouds. 

Our new observations of Q0302 were made to improve both the wavelength calibration 
and resolution of the He + absorption, with enough sensitivity to correlate usefully with the 
HI absorption. This is much more informative than just detecting the mean absorption — 
we can explore the relative contributions of clouds and diffuse gas, as well as measuring 
independently the ionizing radiation field and the helium abundance. A substantial optical 
depth, on the order of the whole effect detected at low resolution (Songaila et al. 1995), is 
expected just from the gas accounted for in the discrete clouds already identified as the HI 
Lyman-a forest. Although this possibility can be modeled theoretically in a statistical way 
for low resolution data (Giroux et al. 1995), our higher resolution He + spectrum allows a 
direct detailed comparison between HI and He + line absorption, and hence a much more 
powerful constraint on models. We find significant He + absorption from HI clouds (with 
optical depth of the order of unity) but also comparable He + absorption even in redshift 
intervals where the best Keck spectrum reveals no detectable HI; thus we directly measure 
absorption attributable separately to both the clouds and the diffuse gas. Even with high 
resolution however our spectrum suggests nonzero flux at all wavelengths, which constrains 
the ionizing background spectrum and leads to an upper limit on the density of diffuse gas. 
Absorption from gas near the quasar, where the incident spectrum is known approximately 
from the direct measurement of the quasar spectrum, allows independent constraints on the 
density and helium abundance of the gas. 
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2. Observations and Reductions 

On 3 separate visits in October, November, and December of 1995, Q0302 was observed 
with the Goddard High Resolution Spectrograph (GHRS). The G140L grating and "Dl" 
detector were chosen especially because of their high UV efficiency, but this setup also 
provides good spectral resolution of 0.6A (about one diode). Q0302 is too faint for a direct 
target acquisition with GHRS, necessitating initial acquisitions with the Faint Object 
Spectrograph, followed by an offset into the GHRS "Large Science Aperture" (LSA). This 
scheme yields a target-positioning uncertainty in the LSA that potentially could translate 
into a ~ 1 A uncertainty in the absolute wavelength calibration. (However, empirically we 
find, for example, an offset of only 0.3 A between the expected and observed wavelengths 
for the strong interstellar C II absorption component at 1334.5 A.) 

The requirement to use the LSA results in two complications related to limiting 
geocoronal/airglow contamination. First, a grating tilt was chosen with coverage of 
1240-1525 A that starts safely redward of the very strong geocoronal Lyman-a line. Second, 
01 1304 A airglow cannot be avoided completely as its wavelength is very near to redshifted 
He + Lyman-a for Q0302. Hence, the most useful science observations may be collected only 
during the spacecraft "nighttime" portions of each HST orbit, and the STScI staff helped 
insure that scheduling occurred in a manner to maximize the availability of such dark time. 

Based on the earlier FOC estimates of the UV flux, we expected a net count rate from 
Q0302 significantly lower than the typical GHRS noise background rate. Hence, GHRS 
science data during the nighttime portions of 16 HST orbits were taken using special 
noise-rejection commanding known as "FLYLIM", described in greater detail below. In 
addition, most of 2 additional HST orbits were devoted to routine "ACCUM" GHRS 
observations of Q0302 during spacecraft nighttime. Finally, there were brief periods at the 
beginnings or ends of each HST target visibility period in which the spacecraft was not in 
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Earth shadow; additional GHRS data were collected in routine ACCUM mode during these 
brief spacecraft "daytime" periods. 

Amongst these various datasets, the FLYLIM/nighttime observations have the greatest 
statistical precision due to both their longer total integration time, plus reduced noise 
background, and hence are the focus of most of our analysis. FLYLIM cut the noise 
background count rate by nearly one-half while rejecting only 7.5% of the integrations. 
However, as described below our GHRS observations reveal a UV flux for Q0302 in excess of 
initial FOC expectations (and even somewhat higher than a recent re-calibration of the FOC 
data by Jakobsen 1996). This higher UV flux interacts with the FLYLIM noise-rejection 
scheme in a nonlinear manner that necessitates more sophisticated reductions than are 
standard in the STScI pipeline routines. The additional 2 orbits of ACCUM/nighttime 
observations are essential to confirming the flux-calibration of the FLYLIM/nighttime 
observations. The ACCUM/daytime observations are strongly contaminated, especially by 
geocoronal emission, and therefore not of much direct use for many of the science issues 
addressed here; nonetheless, they provide useful secondary information on airglow line 
profiles, and on the distribution of background noise events. 



2.1. Calibration of FLYLIM Data 

As the proper flux calibration of the FLYLIM/nighttime observations is not 
automatically handled in the standard STScI pipeline reductions, we discuss here in detail 
the relevant offline flux corrections. The GHRS background noise (mainly cosmic ray 

2 A preliminary analysis using the 01 line profile collected from the daytime portions of 
the orbits suggests that 01 contamination is likely to be small in the nighttime observations 
emphasized here. 
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induced Cerenkov radiation) is non-Poissonian, arriving predominantly in short bursts 
of events, and FLYLIM commanding is intended to detect and reject such noise bursts. 
FLYLIM rejects (onboard the spacecraft) any 0.2s sub-integration in which the total counts 
accumulated in all 500 GHRS diodes equals or exceeds some specified number of counts; 
in this case a threshold of 3 counts was selected. This threshold is substantially larger 
than the count rate expected (per 0.2s) from the QSO itself across all 500 diodes, and thus 
sub-integrations at or exceeding the threshold of 3 counts are presumed to be noise and 
rejected. Typical GHRS observing modes spend of order 1/16 of the actual integration 
time monitoring the noise background and FLYLIM is commanded to reject both on-source 
sub-integrations (source+background) and background monitoring (background only) 
sub-integrations with the identical threshold. 

There is thus the potential with FLYLIM for subtle background and threshold 
rejection effects not accounted for in the standard STScI pipeline flux calibration. It is 
helpful to conceptualize the needed offline flux corrections for FLYLIM data as due to two 
principal effects, although of course these two are coupled (and hence require a coupled 
calibration / correction) . 

First, even if there were no background noise, there would be rare (but expected) 
positive counting fluctuations from the source (QSO) that would occasionally exceed the 
FLYLIM threshold, hence causing the rejection of an entire 0.2s sub-integration. Even 
though rare, such large rejected positive source fluctuations will not be accounted for in the 
standard pipeline estimation of the source count rate, and this artificially suppresses the 
standard pipeline inference about the QSO flux. This effect, of course, is wavelength (or 
diode) dependent, as the QSO count rate itself is wavelength dependent. This effect will be 
especially important in those regions of the spectrum in which the count-rate is highest, 
as it is in such wavelength/diode regimes where the largest positive source fluctuations 
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are most likely. (In the case of Q0302, this will be longward of the He + break). The flux 
correction needed to the pipeline-reduced spectrum for this effect is largely a multiplicative 
one; that is, the actual flux at each wavelength is (approximately) some fixed factor higher 
than that implied by the standard pipeline reductions, and may be simply estimated as 
described below. 

A second, and perhaps more subtle effect, is that the background count rate may also be 
misinterpreted with FLYLIM. For concreteness, consider a case in which source+background 
are being observed, and a background noise event generates 2 counts in a particular 0.2s 
sub-integration while the source generates 1 count in that same sub-integration. This 
0.2s sub-integration on source+background will be rejected by FLYLIM as it attains the 
threshold value of 3. However, when monitoring only the background (no source), there 
would not have been a FLYLIM rejection as the 2 count background event itself falls below 
the rejection threshold. Thus, comparatively lower background rates may lead to FLYLIM 
rejections when observing source+background than those that lead to FLYLIM rejections 
when monitoring background only. The result is that for observations taken in FLYLIM 
mode, the noise background count rate estimated by the standard pipeline reductions (from 
background-only monitoring data) appears higher than the noise background count rate 
that is actually allowed through when observing both background+source. This latter effect 
essentially necessitates a zero-point count-rate (or flux) offset correction to the pipeline 
reduced data. This zero-point offset is especially important for the low count rate portions 
of the FLYLIM spectrum (shortward of the He + break for Q0302), and is corrected for as 
described below. 



2.2. Empirical FLYLIM Flux Correction 



With the above simple conceptual model in mind, an empirical flux correction can 
be made to the standard pipeline reductions of the FLYLIM/nighttime observations, 
by requiring that the corrected spectrum have both the same overall (low-resolution) 
spectral energy distribution and the same mean flux as measured directly from the 
routine ACCUM/nighttime observations. In practice, this is done iteratively. First, an 
initial estimate of the mean zero-point offset flux correction (corresponding to the second 
effect discussed in section 2.1) is obtained from the difference between mean fluxes of 
the routine ACCUM/nighttime and FLYLIM/nighttime pipeline-reduced spectra; the 
mean fluxes compared here are those shortward of the He + break in Q0302, where this 
zero-point offset dominates. Then, with this initial zero-point correction applied, the 
ratio of mean fluxes between the (now partially corrected) FLYLIM/nighttime spectrum 
and the routine ACCUM/nighttime spectrum is determined to estimate the (mainly) 
multiplicative correction corresponding to the first effect discussed in section 2.1; the mean 
fluxes compared here are those longward of the He + break, where the first effect dominates. 
This procedure is iterated until the difference in mean fluxes below the break converges 
to near zero (to within ~ 3 x 10~ 18 erg/sec/cm 2 /A) and simultaneously the ratio of mean 
fluxes longwards of the break converges to unity (to within ~ 0.3%). [Note that for these 
data, in practice it makes no significant difference whether such an empirical correction is 
done in flux or count/rate units, although the latter is more rigorously consistent.] 

The empirically corrected FLYLIM/nighttime spectrum is displayed in figure 1, with 
the accompanying formal error spectrum. The empirical correction in this case may be 
thought of as roughly equivalent to a zero-point offset correction of about 2.5 xl0~ 17 
erg/sec/cm 2 /A coupled with a 30% multiplicative correction to the flux. The error 
spectrum displayed in figure 1 is appropriate for assessing the statistical significance of 
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spectral features, but our empirical correction to bring the FLYLIM/nighttime spectrum 
into agreement with ACCUM/nighttime spectrum also must reflect the zero-point flux 
uncertainty of the ACCUM/nighttime data. Hence it should also be emphasized that the 
empirically corrected FLYLIM/nighttime spectrum of figure 1 has inherited a zero-point 
uncertainty of about 1.1 xlCT 17 erg/s/cm 2 /A (and this zero-point uncertainty is not 
reflected in the formal error spectrum). 

On the other hand, it should also be realized that for most of those discussions below 
that refer to flux measures-e.g., estimates of the optical depth over broad bands due to 
He + absorption-one can use either the entirely routine ACCUM/nighttime spectrum or 
the empirically corrected FLYLIM/nighttime spectrum of figure 1. Both spectra yield very 
similar results (including errors, which are dominated over broad bands by the zero-point 
uncertainty associated with the ACCUM/nighttime observations). 



2.3. Model FLYLIM Flux Correction 

We have also independently confirmed the FLYLIM/nighttime flux corrections by 
a more detailed quantitative model of the FLYLIM rejection process; in contrast to the 
empirical correction of section 2.2, the model correction described now does not depend 
on any recalibration to the ACCUM/nighttime flux. For this more detailed confirmatory 
model, we assume that the source (QSO) count rate (per 0.2s sub-integration) may be 
modeled as a Poisson process. For the distribution of GHRS background noise events, 
a Poissonian model is known to be inadequate so we instead estimate this noise count 
distribution (i.e., what fraction of noise bursts yield 0,1,2,3,... count events in each 0.2s 
sub-integration) from the background-only monitoring data collected during each HST 
science orbit in routine ACCUM mode (i.e., no source and no FLYLIM rejection). These 
background-only monitoring data are essentially background "spectra" in which all 500 
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GHRS diodes are read-down in exposures lasting about 7 sec each. We rebin each such 
background spectrum into bins of ~14 diodes in width, and count the number of background 
events yielding 0,1,2,3,... counts in each bin. Effectively, each bin then samples 14-diodes 
x 7 sec m 100 diode-secs worth of background monitoring data; for the model, we assume 
that this binning scheme yields a similar distribution of background events as that actually 
sampled (aboard the spacecraft) in each 0.2s sub-integration across all 500 diodes, as it is 
also the case that 500-diodes x 0.2s ~ 100 diode-secs. Finally, we also require that the 
model account for the empirical constraints that 7.5% of all integrations with FLYLIM 
implemented were rejected (with the chosen threshold of 3), and that the pipeline-inferred 
(i.e., uncorrected) background count rate while monitoring background-only (no source) 
with FLYLIM is 0.0036 cnts/s/diode. 

The pipeline-reduced FLYLIM/nighttime spectrum is input to the model as an initial 
estimate of the actual spectrum, and the model provides a first estimate of the required flux 
correction (accounting for both the coupled FLYLIM effects described in section 2.1). This 
partially corrected spectrum is then input again into the model to derive an improved flux 
correction, and the procedure iterated until convergence. It should be emphasized that the 
ACCUM/nighttime flux does not in any fashion enter into this model for the flux correction 
to FLYLIM data. 

This model-corrected FLYLIM/nighttime spectrum agrees very well with the 
empirically-corrected FLYLIM/nighttime spectrum of figure 1. The difference between 
the mean fluxes for model- and empirically-corrected FLYLIM/nighttime spectra is only 
l.OxlO -17 erg/s/cm 2 /A, and the difference in means shortward of the He + break is even 
smaller at 6.0xl0 -18 erg/s/cm 2 /A. In summary, the excellent flux agreement between 
empirically-corrected and model-corrected FLYLIM/nighttime spectra on the one hand, 
plus the agreement between the model-corrected FLYLIM/nighttime spectrum and the 
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entirely routine (and empirical) ACCUM/nighttime spectrum on the other hand provide 
strong confirmation that the somewhat subtle offline flux corrections described in this 
section are properly accounted for in the spectrum of figure 1. 

3. Results and Interpretation 

3.1. Summary of Absorption Features 

The spectrum clearly confirms the main result of Jakobsen et al's FOC study: the 
existence of an absorption edge of significant optical depth. It also appears that the 
absorption is indeed mainly produced by He + 304 A line absorption and not by an unrelated 
\ow-z HI Lyman limit cloud. First, this is confirmed by the location of the break, which in 
the GHRS data occurs within about 5.3 A shortward of the quasar systemic redshift (this 
may be compared with the ±20A agreement uncertainty window allowed by the original 
FOC data). Scaling the arguments presented in Jakobsen et al. (1994), the likelihood of 
a chance superposition of an HI Lyman limit system (with column exceeding 10 18 cm -2 ) 
occurring this close shortward of the QSO systemic redshift is only about 0.3%. Second, 
in any case our GHRS spectrum does not reveal HI Lyman series absorption, which would 
have been detectable if a Lyman limit system were responsible for the break. 

Most significantly, we find significant coincidence in redshift of helium and hydrogen 
absorption features — the He + Lyman-a forest as well as the Gunn- Peterson effect. The 
portion of the new spectrum of greatest interest is displayed in figures 2a and 2b, overlaid 
with synthetic spectra, which are models of the absorption expected from the clouds seen in 
the HI Lyman-a forest absorption. Starting with a Keck spectrum of the quasar (courtesy 
of A. Songaila, E. Hu and L. L. Cowie), we have fitted Voigt profiles (using VPFIT, 
Webb 1987) to all the HI line features down to a threshold column density of 10 12 cm -2 
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(e.g., a central optical depth of 0.05 for b = 15km/s), then used the column densities, 
Doppler parameters and redshifts of these clouds to predict the He + absorption, degraded 
to the resolution of GHRS. (The reason for this roundabout modeling, rather than simply 
multiplying the observed HI transmission by a constant factor [as in e.g. Songaila et al 
1995] is to allow us to model either turbulent or thermal cases, depending for each argument 
on whether maximum or minimum helium absorption is appropriate). 

Note that several He + absorption features correlate in detail with features predicted 
from HI; the main edge is itself one example of such a feature, caused by an identifiable 
complex of HI clouds. This edge is centered at z — 3.268, which is less than the systemic 
redshift of the quasar but matches that of the clouds. This is expected since the ionizing 
radiation from the quasar reduces absorption from the diffuse foreground gas more than 
from saturated clouds (Zheng and Davidsen 1995, Giroux et al 1995) The traditional 
redshift z = 3.285 of Q0302, which would place the edge at 1302A instead of the observed 
1296. 4A is probably quite close to the "true" cosmic or systemic redshift of the quasar, 
as nearly the same redshift is derived from both Lyman-a and from other, narrower lines 
including the [OIII] 5007 emission line; see Espey et al. 1996. The coincidence (to within 
~ 1 resolution element) of the three reddest absorption features with the HI predictions- 
especially the main edge- suggests that the absorption indeed arises from high redshift 
helium even though the edge is not at precisely the quasar systemic redshift. 

Beyond the main edge, a significant "shelf" of flux (mean optical depth about 0.8, level 
11.8 ± 1.4 x 10~ 17 erg/sec/cm 2 /A) extends for about another 10A, to about 1283A (20A 
shortwards of the "expected" edge). We attribute this entire shelf to the fact that in this 
region more of the He + is doubly ionized due to the hard radiation from the quasar (the 
"proximity effect"), greatly reducing absorption from diffuse gas. 

The second edge at 1283A coincides with the highest HI column cloud measured, 
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log N(H I) = 15.8, which may even be optically thick in the He + Lyman limit. At shorter 
wavelengths the flux decreases further beyond the second edge. For broad comparison 
with the results of other studies, we note that here the total He + optical depth (which 
includes both cloud and diffuse gas contributions) is T tota i ~ 2.0, if we compare the mean 
flux averaged over the 40A region shortward of the "shelf to a similar region longward of 
the He + break; the 95% confidence range on r tota i is 1.5 to 3.0, compared with the shelf 
average of 0.8. Since the contribution from highly saturated clouds is little affected by the 
proximity effect, the change or extra absorption, an optical depth of the order of unity, 
is probably caused by unsaturated, redshift-space-filling, diffuse intergalactic gas. In this 
sense we confirm quantitatively the conclusions of Jakobsen et al. (1994). 

There is some flux detected even near the very bluest portions of the spectrum of 
figure 1. For example, in a 20 A bin between 1240 and 1260 A the mean flux is 5.0x10 
erg/s/cm 2 /A with a formal (photon counting) error of ±0.6 xl0~ 17 erg/s/cm 2 /A (and the 
zero-point offset error of l.lxl0~ 17 erg/s/cm 2 /A discussed in section 2.2). This is at least 
indicative of a flux recovery towards lower redshift: for example, the mean flux in the 
1240-1260 A bin is marginally larger by about 2.9a than that (2.8 xlO -17 erg/s/cm 2 /A 
with similar errors) in a redder adjacent bin extending over 1260-1280 A. Note that the 
zero-point offset error is not relevant when comparing the difference in flux between these 
adjacent bins. Such a flux recovery is expected in many models and is certainly expected 
from the results of Davidsen et al. 

The predicted spectra in figure 2 depend on the adopted column densities and Doppler 
parameters of the He + . The models displayed in figure 2a simply assume two constant 
ratios for all the clouds, rj = N(He + )/N(HI) = 20 (corresponding to a spectral slope 
a = 1.8, the typical value observed in low redshift quasars and probably appropriate for the 
unabsorbed ionizing spectrum in the near-quasar zone here), and 77 = 100 (corresponding 



-16- 



to our limit a > 3 derived below for the slope far from the quasar). Figure 2b shows the 
result for maximal absorption from an extremely soft spectrum, in order to demonstrate 
that the observed HI clouds alone cannot for any spectrum explain the absorption in the 
Lyman-o; voids. All the models shown in figure 2 have £ = huel^H = 1 for all components 
(pure turbulent broadening, giving the maximal helium/hydrogen optical depth). The 
absorption is quite sensitive to £ but not to rj as the helium absorption by clouds is highly 
saturated; even very large values of rj do not yield significantly more absorption than the 
models shown if £ = 0.5. Due to variations in £ a precise agreement of the simple model 
with the data is not expected even where the spectrum is uniform. 

In the following subsections we will develop several separate quantitative lines of 
argument based on different features of the spectrum. In the "shelf" region of the spectrum 
the absorption is, within our errors (including the uncertainty in the intrinsic quasar 
spectrum), accounted for entirely by the gas in clouds, provided bn e ~ b H for at least 
some of the clouds. Diffuse absorption of large optical depth in this region is however 
not consistent with the data, and this fact can be used together with the known ionizing 
spectrum of the quasar to set an upper limit on the diffuse gas density near the quasar. 
Knowledge of the spectrum yields broad constraints on the absolute value of the helium 
abundance. Shortwards of 1283 A, the forest clouds are not sufficient to explain all of the 
absorption; indeed there is significant helium opacity even in the well-known "void" in the 
HI forest (Dobrzycki & Bechtold 1991) from which we derive a lower limit on rj and and a 
constraint on the intergalactic ionizing spectrum far from the quasar, which roughly agrees 
with the absolute He + ionizing flux estimated from the proximity effect. In spite of the soft 
spectrum there is residual flux throughout the spectrum, which gives a separate upper limit 
on the diffuse gas density far from the quasar and also indicates that it is likely that the 
intergalactic helium is mostly doubly ionized by this epoch. 
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3.2. Conditions Near the Quasar 

In the region near the quasar, the He + -ionizing flux can be estimated directly from 
the observed quasar flux, with several caveats. If there is dust absorption along the line of 
sight, the true continuum level is everywhere greater than observed, and the slope of the 
spectrum is shallower; in addition, the accumulated Lyman limit absorption of foreground 
HI clouds reduces the observed flux shortwards of 912 A in the rest frame (e.g., Vogel 
and Reimers 1995); and the quasar flux may vary on a timescale long compared with 
observational baseline but short compared to relevant ionization and recombination rates. 
These effects lead to uncalibrated uncertainties which can only be removed by additional 
observations; variability and absorption alter the estimate of the local ionizing flux, leading 
to changes in the proximity arguments, while the unknown intrinsic continuum slope leads 
to an uncertainty the continuum flux and hence in all our estimated optical depths. Our 
approach here is to "parametrize our ignorance" in correction factors for these effects. 
The factors can be constrained using other data, for example spectrophotometry of the 
quasar over a larger wavelength range (extending longwards of 912 A in the rest frame), or 
statistical surveys of other objects. 

We assume that the flux at the He + ionizing continuum edge Vi (\ =228 A rest) would 
if unabsorbed be close to the observed flux longwards of 304 A. (The modest equivalent 
width inferred for the 304 A emission line seems to indicate that the helium-ionizing 
flux from Q0302 is largely escaping, and in the 228-304A region the spectrum is likely 
to be be fairly flat.) We allow for the possibility of other sources of absorption along 
the line of sight, such as an accumulation of Lyman continuum absorption, reducing 
the flux by a total factor i? _1 ; the flux at the continuum edge if unabsorbed is then 
fx ~ 2.6 x 10~ 16 i? erg cm -2 s _1 A _1 , (or /„ 8 x 10~ 30 i? erg cm" 2 s~ 1 Hz~ 1 ). We assume 
an intrinsic power law spectral energy distribution f u = f\cv~ 2 cx v~ a . Zheng and Davidsen 
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(1995) estimate that the intrinsic a ps 1.5, although a ps 2.4 has been measured directly 
and a w 1.8 appears to be typical of nearby samples. For some arguments it is appropriate 
to normalize to the observed flux at HI Lyman-a, in which case it is appropriate to use a 
different quantity the multiplicative difference in the foreground absorption between the HI 
and He + Lyman-a, -R304/1216 ~ 4 2 ' 5_a . For the plausible range of a, -R304/1216 is likely in the 
range between 1 and 4. Similarly for some arguments it is appropriate to use the ratio at 
the continuum edges, -R228/912- 

At a point along the line of sight at redshift zq — Sz, the quasar spectral density is 
(Bajtlick et al. 1988) 

F v = Af v Sz- 2 {l + z) 5 [(l + z)^ 2 -l] 2 

if n — 1, and 

F u = U5z' 2 z 2 [l + {z/2)] l l 2 {l + zf 

if fl = 0, independent of H ; for the present case zq = 3.285, we find (F v / f u )Sz 2 = 6.5 x 10 3 
and 6.3xl0 3 respectively. Expressing the offset Sz from the quasar systemic He + Lyman-a 
redshift in terms of observed wavelength offset 5 A = 304^5^ from the quasar He + Lyman-a, 
304A (1 + zq) = 1302A , we estimate the spectral flux 

F v = 1.2 x 10- 23 (5\/20A)- 2 R erg cm" 2 s _1 Hz _1 , 

the ionizing photon flux 

roo 

F y = dv(F v /hu) ~ 1.9 x 10 4 a" 1 (5A/20A)- 2 J Rcm- 2 sec" 1 
and the He + ionization rate (e.g., Osterbrock 1989) 

POO 

T He+ = / dv(F u a u /hv) ~ 3.0 x 10~ 14 (a + 3)^ (5 A/ '20 A)~ 2 Usee' 1 , 
where the photoelectric cross section of He + , a u = 1.6 x 10 _18 (z/j/z/) 3 cm 2 . 
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3.3. Quasar Lifetime and Ionizing Background from the Proximity Effect 

A quasar can doubly ionize all the helium out to a distance of r = [5z/(l + z)][c/H(z)] 
if the photon flux F 7 is maintained for a characteristic time £q = rn# e /3-F 7 , where for 
cosmic abundance (0.08 by number) the density of helium atoms is related to the density 
parameter of gas by rise = 9 x 10~ 7 Q g h 2 (l + z) 3 cm~ 3 . f\ For Q = 1 we estimate a time for 
double-ionizing to offset S A, 

t Q > 0.6 x lO 7 /^ 1 yr aiT 1 (6 A/20 A) 3 {tt g h 2 /10~ 2 ), 

and about twice this for Q = 0. It is therefore plausible for the helium ionization zone to 
extend about 20A shortwards of the quasar redshift if the lifetime of this quasar is of the 
order of the characteristic Eddington/Salpeter evolution time for accretion, 4 x 10 8 eyr for 
typical radiative efficiency e = E tota i/Mc 2 ~ 10%, even if the gas begins as mostly singly 
ionized. 

If we take 5X ~20A as the point where the quasar ionizing flux is equal 
to the ionizing background, the ionizing spectrum has a specific intensity 
J228 ~ 10~ 24 i? erg cm -2 s _1 Hz _1 sr _1 , implying a soft spectrum, with ratio of hydrogen to 
helium intergalactic ionizing fluxes S = J912/J228 ~ 10 3 -R _1 «/9i2,-2i- 

The predictions based on observed quasar populations are quite sensitive to the 
assumed typical value of a and to models of the absorption. Haardt and Madau (1996) 
predict rj rj 40 at this redshift, although the prediction is likely to increase as a is revised 
(from 1.5 say to 1.8; P. Madau, private communication). The flux here corresponds to 

3 We ignore helium recombinations for this argument. In gas with density contrast p/p, 
the recombination rate is n e a>H e + ~ 8 x 10 -17 f2 s /i 2 (p/p)sec -1 which is slower than the Hubble 
rate H fa 3 x 10~ 17 /isec -1 for Q g w 0.01. If much of the gas is in density concentrations, 
recombinations can be important, requiring a longer quasar lifetime. 



-20- 



r] fx 1.7S fx 200 for R = 4 and Jgu-21 = 0.5; within the uncertainties of both arguments, 
we consider the agreement satisfactory, almost remarkable since it is derived from just 
one quasar. (Note the independent limit derived below, r\ > 100, without reference to the 
proximity effect.) 

3.4. Diffuse Gas Near the Quasar 

The lack of diffuse absorption leads to an upper limit on the diffuse gas density in the 
He ++ bubble, since the ionizing flux from the quasar is known. Continuous absorption 
of optical depth t gp is produced by diffuse He + (which in this context means atoms 
producing unsaturated absorption at redshift z), where the density is given by the standard 
Gunn- Peterson formula (Peebles 1993) 

n He + = (87r/3)r GP [A a (l + z a )\^H{z)h.-\ 

where the transition rate for helium Lyman-a, A a (He) = 16A a (H) = 1.0 x 10 10 sec^ 1 . At 
z = 3.285 this becomes n He + = 0.9 x 10^ 9 hr G pcm~ 3 or equivalently 

Q g h 2 = 1.7 x 10- 5 hT GP (n He+ /n H e)~ 1 - 

The fraction of helium in He + is given by 

n H e+/n He = n e a He /T He+ fx 3.6 x 10' 2 {p/ p)Q g h 2 {5\/20A) 2 {a + ^RT 1 ^ 1 ' 2 , 

where the density contrast of the material is (p/ p) and the recombination coefficient is given 
near T 4 = T/WK = 1 by (Spitzer 1978) 

a He + = 3.4 x 10" 13 Z 2 T 4 " 1/2 cm 3 sec" 1 . 

This results in a limit on the diffuse gas density in the bulk of velocity space (and hence, 
the bulk of the intergalactic spatial volume), 

n g h 2 = 2.1 x 10- 2 T^{p/p)-^ 2 (8X/20A)-\a + Zy^R^Tl^h 1 / 2 
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where tqp is the limit on the continuous opacity at the offset 5X. We conservatively estimate 
from our data r GP < 1.35 at « 1285A (a 95% confidence limit, conservatively neglecting 
any discrete cloud opacity contribution), leading to an upper limit on diffuse intergalactic 
gas of 

n g < 0.019i? 1 /2( Q; /l.5)-V2( /l / J)l/2 ) 

tied to the observed quasar flux. 

This limit is interesting since it is of the same order as the baryon density required by 
Standard Big Bang Nucleosynthesis (Walker et al 1991, Smith et al 1993, Copi et al 1995, 
Sarkar 1996, Hogan 1997). Current estimates of the total range from il b & 0.01(/i/0.7)~ 2 
(e.g. Rugers & Hogan 1996) to tt b « 0.05(/i/0.7)- 2 (e.g. Tytler et al. 1996). Realistic CDM 
models predict that most of the baryons should be concentrated in clouds by this epoch, so 
our result accords with expectations even for large baryon density (Croft et al. 1997). 

3.5. Helium Abundance 

Since both hydrogen and helium are mostly ionized, a precise abundance measurement 
is not possible from absorption which only studies a small fraction of the material. However, 
because of the unique information on helium at high redshift (before the bulk of baryonic 
material had even formed into stars), it is interesting to ask how our data quantitatively 
constrains the abundance, even if imprecisely. 

An absolute helium abundance can be estimated from HI and He + Lyman-a absorption 
but only if: (1) the helium abundance is uniform; (2) He + and HI are in ionization 
equilibrium; (3) we know the shape of the ionizing spectrum; (4) absorption is unsaturated, 
so column densities of absorbing species can measured. The total columns of He + and HI 
are then both proportional to the same line integral / dln 2 e with coefficients depending on 
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the abundance and the ionizing spectrum. (Note that even this statement applies only 
to the redshift integrated column densities; there is only agreement at particular redshift 
in the case of negligible thermal contributions to the atomic velocities. Note also that in 
principle, another test is possible: in regions where the ionizing spectrum is uniform, even 
if it is not known, constraints on variations in rj translate into constraints on the spatial 
variations of Y .) 

These conditions are certainly not met here in detail, but the data also do not allow 
arbitrary variations in abundance. For example, ionization equilibrium relates the value of 
the (redshift-) integrated column density ratio r\ to the absolute abundance of helium, 

Y v 

0.24 1-7S' 2 28/912 

provided we know the ratio of ionizing fluxes 5*228/912 • We can explore the range of allowed rj 
by our family of models based on the HI absorption, and the value of 6228/912 * s constrained 
in the region dominated by quasar radiation, subject to the uncertainties discussed above. 
Using the above estimates of quasar fluxes, 5*228/912 ~ 28i? 22 1 8 y 912 , yielding 

-21 « 0.84 (-1) ( R ^Y 

0.24 uo; v 4 j 

The best guesses from the current data are that 10 < rj < 20 and 2 < -R228/912 < 4, so 
that Y must lie within a factor of a few of the standard big bang prediction. Although the 
constraints on rj and -R228/912 will both improve (the first from better signal to noise, the 
second from better spectrophotometry), it is still unlikely that a reliable estimate can be 
made much more precise than this because of the many assumptions required. The most 
interesting new result here is the approximate concordance with Big Bang predictions at a 
large distance and an early epoch, and over a large volume of space. 
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3.6. Diffuse Gas in the HI Lyman-a "Voids": Lower Limit 

Assuming He + is the dominant species leads to a conservative lower limit on the density 
in absorbing gas, 

Vt g h 2 = 1.65 x 10" 5 /ir GP (n(iJe + )/n(iJe))- 1 [(l + 2)/4.285]" 3/2 . 

In the HI Lyman-a void redshift range (especially near X HI m 4xl266A), the optical depth 
of He + absorption required by our data, after allowing for the absorption from identified 
clouds (which in the void is almost independent of 77), is still r GP > 1.3 (95% confidence), 
requiring a diffuse density, 

Q 9 >3x 10' 5 (h/0.7)-\ 

This is essentially the same as derived by Jakobsen et al., except that we can now rule out 
the possibility of producing this opacity with HI clouds down to the detectability threshold 
t~gp < 0.05 in HI absorption N(HI) = 10 12 cm~ 2 . The higher resolution here shows that the 
helium opacity appears even between the most rarefied detected HI clouds. Of course for 
77 > 100, much of the He + absorption could still be from saturated lines not yet resolved in 
our data. 

In current models, absorption is produced in components of lower column density but 
these components are produced by "clouds" which are indistinguishable from (are really just 
parts of) the diffuse, space-filling protogalactic medium (Cen et al. 1994; Hernquist et al. 
1996; Croft et al 1997; Zhang et al 1997; Bi and Davidsen 1997). The opacity required by 
our data is roughly in accord with these models. We will argue below that it is implausible 
to evacuate space with very high efficiency so that the bulk of the helium must be doubly 
ionized, and even then the constraints on the spectrum impose an interesting upper limit 
on the density. 
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3.7. Ionizing Spectrum Far from the Quasar 

The large ratio of helium to hydrogen optical depths in the void indicates a soft 
ionizing spectrum far from the quasar. If the HI and He + are both optically thin, the 
ratio of optical depths can be used to constrain 77 directly via r GP (He + ) /t gp (HI) = rj/i 
(see Miralda-Escude 1993). In ionization equilibrium 77 is related to the spectral softness 
parameter S (e.g., Giroux et al. 1995), 

a H 1 He+ n H 

where a and T are the recombination and ionization rates respectively for the two species, 
and we have assumed the cosmic abundance 0.08 and a temperature T 4 = 2. In the 
Lyman-a void redshift range (i.e., near Xhi ~ 4 x 1266A), the average optical depth of 
diffuse HI absorption allowed is at most 0.05. (Since uniform HI absorption could have 
escaped detection at this level, this limit of 0.05 assumes that there are some variations in 
t(z), as expected from simulations.) The optical depth of He + is at least 1.3, requiring 
i] > 100, S > 63, and hence an ionizing spectrum of a > 3.0. (This limit becomes stronger 
if the helium absorption is not from uniform gas or comes from saturated absorbers.) 
This is consistent with only soft radiation having reached this material as expected if He 4 " 
absorption is still strongly modifying the emitted spectra of quasars. The evidence for a 
soft spectrum is consistent with that inferred from SilV/CIV ratios (Songaila and Cowie 
1996, Savaglio et al 1997), especially considering the sensitivity of these estimates to details 
of the spectrum and relative metallicity (Giroux and Shull 1997). 

We cannot (without a model of the gas distribution) derive from this data an upper 
limit to i] far from the quasar: if we assume only the minimal (thermal) absorption from 
the clouds, even r\ m 5000 does not yield excessive absorption from the clouds. 
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3.8. Upper Limit on z— filling Gas from the HI voids 



Although we know that there is diffuse gas between the forest clouds, we also know 
there cannot be too much of it or else there would be no light getting through, whereas we 
have an upper limit of t gp < 3 on the mean optical depth. Since we also know the ionizing 
spectrum is soft, we can deduce a limit on the IGM density, tied not to the ionizing flux 
from the quasar (as we did above) but to the cosmic ionizing flux at the HI Lyman edge, 
^912,-21 (i n units of 10 _21 erg cm -2 s _1 sr _1 Hz _1 ), which has other observational constraints 
such as the HI clouds proximity effect (eg Madau and Meiksin 1994): 



We thus get a conservative limit Q g = 0.018 by taking two 95% limits, one to constrain 
i] > 100 (from the lower limit r > 1.3 in the void, above) and one from the upper limit on 
the mean total r < 3 everywhere, and using typical estimates Jgi2,-2i ~ 0.5 (Haardt and 
Madau 1996, Giallongo et al 1996, but see Cooke et al 1996). For a better estimate we 
should allow for the absorption we know is coming from the HI clouds; from our near-quasar 
analysis, we guess that r ~ 1 from identified clouds, allowing only 1.3 < r < 2 more from 
^-filling gas. A reasonable guess for the diffuse gas density is then il g 0.01(/i/0.7)~ 1,5 , 
or even less if rj is larger than 100 as suggested by the proximity effect. The limit is 
significantly better than that from studies from the HI Gunn- Peterson effect, which yield 
limits il g < 0.2J 912 -2i at z = 3 (Giallongo et al. 1992), and comparable to the (more 
model-dependent) limit Q g < 0.01 at z — 4.3 (Giallongo et al 1994). 




which yields 



Q g = 0.018(r GP /3)°- 5 (V0.7)- 1 - 5 (Vl00)- a5 (J 912i _ 21 /0.5) 



0.5 
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3.9. Ionization History 

Although it is possible that He + is the dominant species in the intergalactic gas, the 
upper limit of 3 on the mean optical depth would allow in this case at most a density 
of fig < 7 x 1CT 5 in diffuse gas — a number so low that it appears more likely that the 
ionization bubble around Q0302 is a "proximity effect" — that is, double ionized helium is 
already predominant everywhere, and we are just seeing the region nearest the quasar with 
an even higher ionization. The main features of our observed spectrum are indeed predicted 
by models of the proximity effect in which a significant contribution to the mean He 4 " 
opacity comes from the forest clouds (Giroux et al. 1995). [Note that HS 1700+64 shows no 
such proximity effect (Davidsen et al. 1996), which is most easily explained if the absorption 
somewhat later at z — 2.72 is everywhere dominated by saturated lines in clouds.] The 
comoving radius today of the observed He ++ bubble is H r = [(5X/20A)] x 4600km/sec 
and 2200 km/sec respectively for open and flat cosmologies. If the Q0302 bubble is typical, 
the protogalactic gas has entropy and ionization-state inhomogeneities of the order of 
unity on this scale. The scale of the bubble is not negligible compared to the scale over 
which the "recovery" appears to occur in our spectrum, so it is not clear whether we are 
seeing a cosmic trend or merely the history of radiation percolation (to the redshift of the 
"recovery") along this line of sight. There is thus strong motivation for obtaining a high 
quality spectrum of Q0302 that extends to lower redshifts than probed by our GHRS data. 

4. Summary 

The absorption observed here is broadly consistent with the expectations of hierarchical 
models of structure formation and with conservative models of cosmic ionization. It is clear 
that this type of data will be an important constraint on models and their parameters, 
especially concerning the most diffuse gas filling the bulk of spatial volume. 
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The main new conclusions from the current data are: 1. The He + Lyman-a forest is 
detected; 2. The "diffuse" (redshift-space-filling) medium is also detected, and must have a 
low density (Q < 0.01(/i/0.7)~ 3 / 2 ) consistent with standard primordial nucleosynthesis and 
models of early gas collapse into protogalaxies; 3. The intergalactic ionizing spectrum is soft 
(rj > 100), although the intergalactic helium is probably mostly doubly ionized by z = 3.3; 
4. The helium abundance is within a factor of a few of standard Big Bang predictions, over 
a large volume of space at high redshift. 

There is clearly a strong motivation to get a spectrum of other quasars of comparable 
quality at the same redshift; if we are to draw universal generic conclusions about cosmic 
ionization history, it would be prudent both to extend our results to lower redshift in 
Q0302, and to have more than a single line of sight to check assumptions about intrinsic 
quasar properties and uniformity on different sightlines. In spite of the persuasive checks of 
the GHRS calibration, it would also be good to verify that the zero level is correct in order 
to strengthen our limit on the diffuse gas density. These programs are now underway with 
HST/STIS. 



5. Acknowledgments 

We are grateful for advice and encouragement from A. Davidsen, P. Jakobsen and B. 
Margon, for technical support from A. Suchkov, A. Berman, C. Leitherer and others at 
STScI, for the Keck spectrum from A. Songaila, E. Hu and L. L. Cowie, and for useful 
comments from B. Espey and S. Heap. This work was supported at the University of 
Washington by NASA, and is based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is operated by AURA, 
Inc. under NASA contract. 



-28- 

REFERENCES 

Bajtlick, S., Duncan, R. C. and Ostriker, J. P., 1988, Ap J 327, 570. 
Bi, H. and Davidsen, A. F. 1996, ApJ, in press 

Cen, R., Miralda-Escude, J., Ostriker, J. P., and Rauch, M. 1994, ApJ 437, L9. 

Cooke, A. J., Espey, B. and Carswell, B. 1996, MNRAS, in press. 

Copi, C. J., Schramm, D. N., and Turner, M. S. 1995, Science 267, 192. 

Croft, R. A. C, Weinberg, D. H., Katz, N. and Hernquist, L. 1997, ApJ, submitted 
PtrFph 96110531 ) 

Davidsen, A. F., Kriss, G. A. and Zheng, W., 1996 Nature 380, 47. 

Dobrzycki, A. and Bechtold, J., 1991, ApJ 377, L69. 

Espey, B. and Armus, L., 1996, manuscript in preparation. 

Giallongo, E., Cristiani, S., Trevese, D., et al. 1992, ApJ, 398, L9. 

Giallongo, E., et al. 1994, ApJ, 425, LI. 



Giallongo, E., Cristiani, S., D'Odorico, S., et al. 1996, ApJ, to appear (|astro-ph 9602026 
Giroux, M. L., Fardal, M. A., and Shull, J. M., 1995, ApJ 451, 477. 
Giroux, M. L. and Shull, J. M., 1997, Astron. J., submitted. 
Haardt, F., and Madau, P. 1996, ApJ 461, 20 

Hernquist, L, Katz, N., Weinberg, D. H., and Miralda-Escude, J. 1996, ApJ 457, L51. 
Hogan, C.J., in Critical Dialogues in Cosmology, ed. N. Turok (Princeton), in press 



-29- 

Jakobsen, P., et al., 1994, Nature 370, 35. 

Jakobsen, P., 1996, Science with the Hubble Space Telescope-II, (eds. P. Benvenuti, F. D. 
Macchetto, & E. J. Schreier), STScI, 153. 

Madau, P. and Meiksin, A., 1994, Ap J 433, L53. 

Miralda-Escude, J., 1993, Mon. Not. R. Astron. Soc. 262, 273. 

Miralda-Escude, J., and Rees, M. J., 1994, Mon. Not. R. Astron. Soc. 266, 343-352. 

Miralda-Escude, J., Cen, R., Ostriker, J. P., and Rauch, M., Ap J 471, 582. 

Osterbrock, D. E., 1989, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei, 
University Science Books. 

Peebles, P. J. E., 1993, Principles of Physical Cosmology, Princeton, p. 550. 

Rauch, M., et al. 1997, ApJ, submitted ( |astro-ph 9612245Q 



Rugers, M. and Hogan, C. J., 1996, ApJ 459, LI. 
Sarkar, S., 1996, Rep. Prog. Phys, submitted 

Savaglio, S., Cristiani, S., D'Odorico, S., Fontana, A. Giallongo, E. and Molaro, P. 1997, 
Astron. Astrophys., in press ( |astro-ph 9606063|) 

Smith, M. S., Kawano, L. H., and Malaney, R. A. 1993, ApJ 85, 219 

Songaila, A. and Cowie, L. L., 1996, Astron. J., 112, 335. 

Songaila, A., Hu, E. M., and Cowie, L. L., 1995, Nature 375, 124. 

Spitzer, L., 1978, Physical Processes in the Interstellar Medium, Wiley. 

Tytler, D., Fan, X.-M. and Buries, S., 1996, Nature 381, 207. 



-30- 



Tytler, D. and Jakobsen, P., 1996, manuscript in preparation. 

Tytler, D., et al., 1995, in QSO Absorption Lines (eds. Bergeron, J., Meylan, G., and 
Wampler, J.), 289. 

Vogel, S. and Reimers, D. 1995, Astron. Astrophys. 294, 377 

Walker T. P., Steigman, G., Schramm, D. N., Olive, K. A. and Kang, H.-S. 1991, ApJ 376, 
51 

Webb, J., Ph.D., 1987, thesis, University of Cambridge. 

Zhang, Y., Anninos, P., Norman, M., and Meiksin, A. 1997, ApJ, submitted (|astro-ph 



9609194| ) 



Zheng, W. and Davidsen, A., 1995, ApJ 440, L53. 



This manuscript was prepared with the AAS L^TjrjX macros v4.0. 



-31 - 



Fig. 1. — HST/GHRS spectrum and formal error spectrum (latter appropriate for assessing 
the significance of spectral features), both displayed at the instrumental resolution of 
0.6A, using an iterative empirical flux correction for the nighttime FLYLIM observations. 
The correction scheme leaves a residual la uncertainty in the zero level of 1.1 x 
10~ 17 erg cm -2 s _1 A _1 . 

Fig. 2. — A portion of the HST spectrum overlaid with a model spectrum predicted 
on the basis of the model distribution of HI derived from a Keck spectrum of the HI 
Lyman-a forest. The quasar emission spectrum is fitted with a flat continuum with flux 
2.6 x 10~ 16 erg cm~ 2 s -1 A -1 , plus an He+ emission line (centered at 1303A, FWHM 4000 
km/sec, and equivalent width of 5A). Ticks indicate the fitted HI velocity components from 
the Keck spectrum. Doppler parameters and column densities from the fit were used to 
predict the He + absorption spectrum at the GHRS resolution. Two predictions are shown 
in figure 2a (upper panel), dotted and dot-dash curves corresponding to r] = 20 and 100 
respectively, both models assuming pure turbulent broadening, b He + = bur, figure 2b (lower 
panel) shows r] = 500 (and r] = 100 again for comparison). Note: (1) the HST and Keck 
spectra appear to show corresponding absorption features near the He + edge; (2) rj = 20 is 
probably sufficient to explain the absorption features near the quasar entirely with clouds; 
(3) a large (tgp > 1.35) Gunn- Peterson optical depth is only allowed outside the proximity 
of the quasar, below about 1283A; (4) there is significant He + opacity (tgp > 1.3) even at 
the redshift of the conspicuous HI Lyman-a forest void near 1266A; (5) there is significant 
nonzero flux even far from the quasar. 
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